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ABSTRACT 

This is the first paper of a series presenting a Spheroid's Panchromatic Investigation in Dif- 
ferent Environmental Regions (SPIDER). The sample of spheroids consists of 5,080 bright 
(M r < -20) Early-Type galaxies (ETGs), in the redshift range of 0.05 to 0.095, with optical 
(griz) photometry and spectroscopy from SDSS-DR6 and Near-Infrared (YJHK) photometry 
from UKIDSS-LAS (DR4). We describe how homogeneous photometric parameters (galaxy 
colors and structural parameters) are derived using grizYJHK wavebands. We find no sys- 
tematic steepening of the CM relation when probing the baseline from g-r to g-K, implying 
that internal color gradients drive most of the mass-metallicity relation in ETGs. As far as 
structural parameters are concerned we find that the mean effective radius of ETGs smoothly 
decreases, by 30%, from g through K, while no significant dependence on waveband is de- 
tected for the axis ratio, Sersic index, and a4 parameters. Also, velocity dispersions are re- 
measured for all the ETGs using STARLIGHT and compared to those obtained by SDSS. 
The velocity dispersions are re-derived using a combination of simple stellar population mod- 
els as templates, hence accounting for the kinematics of different galaxy stellar components. 
We compare our (2DPHOT) measurements of total magnitude, effective radius, and mean 
surface brightness with those obtained as part of the SDSS pipeline (Photo). Significant dif- 
ferences are found and reported, including comparisons with a third and independent part. A 
full characterization of the sample completeness in all wavebands is presented, establishing 
the limits of application of the characteristic parameters presented here for the analysis of the 
global scaling relations of ETGs. 
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1 INTRODUCTION 

We have witnessed tremendous advances in our ability to measure 
galaxy properties with unprecedented accuracy over the past two 
decades. The greatest leap forward was the advent of the CCD, with 
an order of magnitude increase in sensitivity over photographic film 
and the ability to easily make quantitative measurements. Today, 
surveys like the Sloan Digital Sky Survey (SDSS) and UKIRT In- 
frared Deep Sky Survey (UKIDSS) provide access to high qual- 
ity data covering a large wavelength range, probing different astro- 
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physical aspects of the galaxies. These advances have been imme- 
diately applied to examining the global properties of galaxies. 

Understanding the formation and evolution of galaxies re- 
quires probing them over a long time (redshift) baseline to establish 
the physical processes responsible for their current observed prop- 
erties. It is far easier to measure nearby ETGs as opposed to their 
counterparts at high redshift, and this must be borne in mind when 
comparing samples at opposite distance extremes. Almost as soon 
as high quality data became available for nearby samples of ETGs 
it was recognized that they occupy a 2-dimensional space, the fun- 
damental plane (FP), represented by the quantities radius, velocity 
dispersion, and surface brightness R- a - fi. Brosche (1973) was the 
first to examine ETGs using multi-variate statistical techniques, ap- 
plying Principal Component Analysis (PCA) to the relatively poor 
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data available. Although his results were not fully appreciated at the 
time, it drew the attention of other researchers who further studied 
the implications of the FP (Djorgovski 1987; Dressier et al. 1987). 
Brosche's fundamental contribution was to show that we should 
be looking for sets of data with the smallest number of significant 
principal components when starting from a large number of input 
parameters. In this way we reduce a high-dimensional dataset to 
only those quantities that are likely to have physical meaning. Sev- 
eral studies followed Brosche's: Bujarrabal, Guibert & Balkowski 
(1981), Efstathiou & Fall (1984), Whitmore (1984), and Okamura 
etal. (1984). 

The fundamental plane is a bivariate scaling law between R, 
a, and I, where fi = -2.5 log I, expressed as R ~ <r A I B . In order to 
obtain accurate and meaningful coefficients that can be compared 
to theoretical expectations (such as those from the Virial Theorem, 
which implies A = 2 and B= -1) we need to not only have a ho- 
mogeneous sample of ETGs, but also to understand the selection 
effects in defining the sample and to properly measure the photo- 
metric and spectroscopic quantities involved. Several contributions 
in the past have met most of these requirements (Graham & Colless 
1997; J0rgensen et al. 1996). However, the lack of homogeneous 
data covering a large wavelength baseline while also probing the 
entire range of environments (local galaxy density) impeded fur- 
ther progress. Bernardi et al. (2003a,b), Bernardi et al. (2006), and 
Bernardi et al. ( 2007) were the first to fill this gap and set strong 
constraints on the FP coefficients and their implications, though 
limited to the optical regime (see also Hyde & Bernardi 2009). 
Another important and often overlooked aspect of such studies is 
the impact of different techniques and implementations for mea- 
suring R - a - n and their respective errors, which ultimately will 
be propagated and compared to the distribution of residuals around 
the FP (Prugniel & Simien 1996; Gargiulo et al. 2009). 

This is the first paper of a series presenting the Spheroid's 
Panchromatic Investigation in Different Environmental Regions 
(SPIDER) survey. SPIDER utilizes optical and Near-Infrared 
(NIR) photometry in the grizY J H K wavebands as well as 
spectroscopic data. Spectroscopy and optical photometry are 
taken from SDSS DR6, while the YJHK data are obtained 
from the UKIDSS-LAS DR4. In the present work (Paper I), 
we describe how the sample of ETGs is selected, how the 
photometric and spectroscopic parameters are derived for each 
galaxy, and derive an accurate estimate of the completeness 
of the sample in each band. The grizYJHK galaxy images 
have been homogeneously analyzed using 2DPHOT, an auto- 
matic software designed to obtain both integrated and surface 
photometry of galaxies in wide-field images (La Barbera et 
al. 2008a). We present a detailed comparison of the 2DPHOT 
output quantities (magnitudes and structural parameters) to those 
provided by the SDSS Photo pipeline (Stoughton et al. 2002). 
We have also re-computed central velocity dispersions from the 
SDSS spectra using the software STARLIGHT (Cid Fernandes 
et al. 2005), and compared these new estimates to those from 
SDSS. Velocity dispersions are re-derived using a combination 
of simple stellar population models as templates. This procedure 
minimizes the well known template mismatch problem, accounting 
for the different kinematics of various stellar components in a 
galaxy. All the photometric and spectroscopic measurements 
presented here are made available through an ascii table at 
http ://w w w. lac . inpe .br/bravo/arquivos/S PIDER_data_paperI .ascii . 1 
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The complete SPIDER data-set is also made available, on request, 
through a database structure which allows the user to easily retrieve 
all information by issuing SQL queries. 

In Sec. 2, we describe how the galaxy sample is selected. 
Sec. 3 describes how grizYJHK-band images are analyzed to 
derive integrated photometry and the structural parameters, with 
corresponding uncertainties. Sees. 4 and 5 compare the overall in- 
tegrated and structural properties of ETGs from g through K, de- 
riving color-magnitude relations and presenting the distribution of 
structural parameters in all wavebands. Sec. 6 compares the struc- 
tural parameters derived from 2DPHOT with those from SDSS. In 
Sec. 7, we describe the measurement of central velocity disper- 
sions, comparing them to those from SDSS. The completeness of 
the sample is studied in Sec. 8. A summary is provided in Sec. 9. 

Throughout the paper, we adopt a cosmology with Ho = 
75 km s" 1 Mpc" 1 , fi m = 0.3, and ft A = 0.7. 



2 SAMPLE SELECTION 
2.1 Sample Definition 

The sample of ETGs is selected from SDSS-DR6, following a pro- 
cedure described in La Barbera et al. (2008b) and La Barbera & de 
Carvalho (2009), selecting galaxies in the redshift range of 0.05 to 
0.095, with 0A M r <- 20, where 01 M r is the k-corrected SDSS 
Petrosian magnitude in r-band. The k-correction is estimated using 
the software kcorrect (version 4_1_4; Blanton et al. 2003a, here- 
after BL03), through a restframe r-band filter blue-shifted by a fac- 
tor (1 + zo) (see also Sec. 3.1). As in previous works (e.g. Hogg et 
al. 2004), we adopt zo = 0.1. The lower redshift limit of the sam- 
ple is chosen to minimize the aperture bias (Gomez et al. 2003), 
while the upper redshift limit guarantees (1) a high level of com- 
pleteness (according to Sorrentino et al. 2006) and (2) allows us 
to define a volume-limited sample of bright early-type systems. In 
fact, ETGs follow two different trends in the size-luminosity di- 
agram (Capaccioli, Caon & D'Onofrio 1992; Graham & Guzman 
2003). The separation between the two families of bright and or- 
dinary ellipticals occurs at an absolute B-band magnitude of —19, 
corresponding to the magnitude limit of M r ~ —20 we adopt here. 
At the upper redshift limit of z = 0.095, the magnitude cut of 
—20 also corresponds approximately to the magnitude limit where 
the SDSS spectroscopy is complete (i.e. a Petrosian magnitude of 
m r ~ 17.8). Following Bernardi et al. (2003a), we define ETGs 
using the SDSS spectroscopic parameter eClass, that indicates the 
spectral type of a galaxy on the basis of a principal component 
analysis, and the SDSS photometric parameter fracDev r , which 
measures the fraction of galaxy light that is better fitted by a de 
Vaucouleurs (rather than an exponential) law. In this contribution, 
ETGs are those systems with eClass < and fracDev r > 0.8. 
We select only galaxies with central velocity dispersion, ao, avail- 
able from SDSS-DR6, in the range of 70 and 420 km s~\ and with 
no spectroscopic warning on (i.e. zWarning attribute set to zero). 
These constrains imply retrieving only reliable velocity dispersion 
measurements from SDSS. All the above criteria lead to a sample 
of 39, 993 ETGs. 

The sample of ETGs with optical data is then matched to the 
fourth data release (DR4) of UKIDSS-Large Area Survey (LAS). 
UKIDSS-LAS DR4 provides NIR photometry in the YJHK 
bands over ~ 1000 square degrees on the sky, with significant over- 
lap with SDSS (Lawrence et al. 2007). The YHK-bmd data have 
a pixel scale of 0.4 " /pixel, matching almost exactly the resolu- 
tion of the SDSS frames (0.396 " /pixel). J band observations are 
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Figure 1. Distribution in RA and DEC of the SPIDER sample. Black 
points mark the optical griz data, while red and blue symbols denote the J 
and YJHK data, respectively. . 



carried out with a resolution of 0.4 " /pixel, and then interleaved to 
a subpixel grid. This procedure results into stacked frames with a 
better resolution of 0.2 " /pixel. The YJHK stacked images (mul- 
tiframes) have average depths 2 of 20.2, 19.6, 18.8, and 18.2 mags, 
respectively. The matching of SDSS to UKIDSS data was done 
with the CrossID interface of the WFCAM Science Archive web- 
site 3 . For each ETG in the SDSS sample, we searched for the near- 
est UKIDSS detection within a radius of 1", by considering only 
UKIDSS frames with better quality flag (ppErrBits < 16). The 
matching result was very insensitive to the value of the searching 
radius. In fact, changing it to 0.5" leads to decrease the sample with 
galaxies having YJHK data available by only five objects. The 
number of matched sources is maximum in J band, with 7, 604 
matches, while amounts to 5,698, 6,773, and 6,886 galaxies in Y, 
H, and K bands, respectively. Considering ETGs simultaneously 
matched with two UKIDSS bands, H + K provides the maximum 
number of objects (6, 575). For any possible set of three bands, the 
number of matches varies between 5, 228 (Y + J + K) to 5, 323 
(J + H + K), which is not significantly larger than the number of 
5, 080 ETGs having photometry available in all the YJHK bands. 
For this reason, we have retrieved NIR data for only those galaxies 
with available photometry in either J (7, 604), or H + K (6, 575), 
or YJHK wavebands (5, 080). The completeness in magnitude of 
each sample is characterized in Sec. 8. 

In summary, the SPIDER sample includes 39,993 galaxies 
with available photometry and spectroscopy from SDSS-DR6. Out 
of them, 5,080 galaxies have NIR photometry in DR4 of UKIDSS- 
LAS. The distribution of galaxies with optical and NIR data on the 
sky is illustrated in Fig. 1. 



2 defined by the detection of a point source at 5<r within a 2" aperture. 

3 See http://surveys.roe.ac.uk/wsa/index.html for details. 



Specification of a given family of systems means setting a property 
(or properties) that isolates systems that presumably went through 
similar evolutionary processes. However, when we select as ETGs 
systems with SDSS parameters eClass < and fracDev r > 0.8 
we expect a certain amount of contamination by galaxies exhibiting 
faint structures resembling spiral arms or other non-systemic mor- 
phologies. Although several morphological indicators have been 
proposed from the parameters of the SDSS pipeline (e.g. Strat- 
eva et al. 2001), the eyeball classification is still considered one 
of the most reliable indicator despite its evident subjectivity (Wein- 
mann et al. 2009). We have visually inspected a subsample of 4,000 
randomly chosen galaxies from our sample, classifying them into 
three groups: ETGs (featureless spheroids); face-on LTGs - late 
type galaxies (a bulge surrounded by an obvious disk); and edge- 
on LTGs (a bulge with a prominent disk). This classification is then 
used to evaluate the ability of the different morphological indica- 
tors to distinguish edge-on and face-on galaxies from the bonafide 
ETGs. 

Five SDSS morphological indicators are considered as shown 
in Figs. 2(a-e): (i) The r-band Inverse Concentration Index (ICI r ), 
defined by the ratio of the 50% to the 90%-light Petrosian radii 
(see Shimasaku et al. 2001); (ii) The parameter fracDeV r , cor- 
responding to the fraction of the total fitted model accounted for 
the de Vaucouleurs component. Notice that the total fit is not a 
bulge+disk decomposition but a re-scaled sum of the best-fitting 
de Vaucouleurs and exponential components (Bernardi et al. 2006); 
(iii) The eClass indicator (see Sec. 2.1); (iv) The fractional likeli- 
hood of a de Vaucouleurs model fit (fLDeV r ), defined as: 



fLDeVr 



InLDeVr 



(InLDeVr + lnLexp r + InLstar) 



(1) 



where LDeV, Lexp and Lstar are the probabilities of achiev- 
ing the measured chi-squared for the de Vaucouleurs, exponential 
and PSF fits, respectively; and (v) The projected axis ratio (b/a) r 
of a de Vaucouleurs fit (deVAB r - an SDSS attribute). From Fig. 2 
fLDeV r and (b/a) r are the two better performance indicators dis- 
criminating ETGs from face-on LTGs (see Mailer et al. 2009). 
Based on the visual inspection, we define two cutoff values (see 
panel f in the Figure), one for each morphology indicator, defin- 
ing a region where the contamination rate is ~5%. We notice that 
this contamination rate is ~2.5 times smaller than that for the sam- 
ple of 4, 000. The selected values for the cutoffs are 0.04 and 0.4 
for fLDeVr and (b/a) r , respectively. The same constraints imposed 
to the entire ETG sample, define a sub-sample of 32,650 bonafide 
ETGs. We flag 7,343 objects as lying in the non-bonafide ETGs 
region denned by (b/a) r < 0.4 or fLDeV,. > 0.04, so that we 
can study the impact of contaminants on the global properties of 
bonafide ellipticals. 



3 PHOTOMETRY 

For each ETG, we retrieved the corresponding best-calibrated 
frames from the SDSS archive and the multiframes from the WF- 
CAM Science Archive. In the case of SDSS, only the griz images 
were analyzed, since the signal-to-noise of the u-band data is too 
low to measure reliable structural parameters. The resulting photo- 
metric system, consisting of the grizYJHK wavebands, is shown 
in Fig. 3, where we plot, for each band, the corresponding overall 
transmission curve. Regarding photometric calibration, the zero- 
point of each image was retrieved from the corresponding SDSS or 
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Figure 2. Different morphological indicators are plotted as a function of the central velocity dispersion of the galaxy for a subsample of 4,000 ETGs (panels 
a-e). Panel (f) plots the two optimal indicators and the dashed lines indicate the region where the contamination is ~6% ((b/a) r > 0.4 and fLDeV r < 0.04. 
Face-on LTGs are indicated in black, edge-on LTGs in blue and bonafide ETGs in red, as shown in the upper-right corner of panel (d). 



UKIDSS archives. As a result of the different conventions adopted 
in the two surveys, the SDSS photometry is in the AB photomet- 
ric system 4 , while UKIDSS data are calibrated into the Vega sys- 
tem (Lawrence et al. 2007). 

3.1 Integrated properties 

We have measured both aperture and total galaxy magnitudes with 
different methods, homogeneously for both the optical and NIR 
data. Aperture magnitudes are estimated with S-Extractor (Bertin 
& Arnout 1996). For each galaxy, several apertures are measured, 
spanning the 2 to 250 pixels diameter range. A set of adaptive aper- 
ture magnitudes is also measured. The adaptive apertures have di- 

4 We actually apply small offsets to the griz zero-points in order 
to produce a better match between the SDSS and AB systems (see 
www.sdss.org/dr6/algori thms/fluxcal.html#sdss2ab). 



ameters of dk = k ■ r Kr t , where k — 3, 4, 5, 6 is a multiplicative 
factor, and r Kr i is the Kron radius (Kron 1998) in the i-band, as 
estimated with S-Extractor. The r Kr i is measured in i-band since 
this is approximately in the middle of the SPIDER photometric sys- 
tem (Fig. 3). For k — 3, the median value of dk is ~ 5.8", which 
is more than five times larger than the typical FWHM in all wave- 
bands (see Sec. 3.2). More than 95% of all the ETGs have a value 
of dz larger than ~ 4.8", and a ratio of ds to the seeing FWHM 
value larger than 3. This makes the adaptive aperture magnitudes 
essentially independent of the seeing variation from g through K 
(see Sec. 3.2). Different types of total magnitudes are adopted. For 
the optical wavebands, we retrieve both petrosian and model mag- 
nitudes, m p and m^, from the SDSS archive (see Stoughton 
et al. 2002). For each waveband, the Kron magnitude, micr, is also 
measured independently for all galaxies, within an aperture of three 
times the Kron radius in that band. The Kron magnitudes are ob- 
tained from MAG-AUTO parameter as estimated in S-Extractor. 
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Figure 3. Transmission curves, t(A), for the grizY J H K filters. Each 
curve has been normalized to an area of one. 



Finally, each galaxy has the estimate of total magnitude, my, from 
the corresponding two-dimensional fitting model (see Sec. 3.2). 

To obtain homogeneous measurements from g through K, 
magnitudes are k-corrected and dereddened for galactic extinc- 
tion by re-computing both corrections with the same procedure in 
all wavebands, rather than retrieving them, when available, from 
the SDSS and UKIDSS archives. For each galaxy, the amount 
of extinction is estimated from the reddening maps of Schlegel, 
Finkbeiner, and Davis (1998), applying the correction of Bonifa- 
cio, Monai & Beers (2000) that reduces the color excess value, 
E(B - V), in regions of high extinction (E(B - V) > 0.1). 
This correction is not included in the SDSS database, and only a 
very small fraction of ETGs (< 1%) is found in the high extinction 
regions. We computed k-corrections using the software kcorrect 
(BL03), through restframe filters obtained by blue-shifting the 
throughput curves in Fig. 3 by a factor (1 + 20). For zo — 0, one 
recovers the usual k-correction. For galaxies at redshift z = zo, the 
k-correction is equal to —2.5 log(l + zo), independent of the filter 
and the galaxy spectral type. We have adopted zo = 0.0725 which 
is the median redshift of the ETG sample. According to BL03, 
this choice allows uncertainties on the k-corrections to be mini- 
mized 5 . We have tested how the waveband coverage can affect the 
k-corrections. For the sample of 5, 080 galaxies with available data 
in all grizYJHK bands, we have estimated the k-corrections in 
griz for two cases, where we used (i) all the eight wavebands and 
(ii) only the SDSS bands, k-corrections turned out to be very stable 
with respect to the adopted waveband' s set, with the standard de- 
viation of k-correction differences being smaller than O.Olmag for 
all the griz wavebands. 



5 The value of zq = 0.0725 is smaller than that of 20 = 0.1 adopted for 
sample selection (Sec. 2). The value of zo = 0.1 makes the selection more 
similar to that performed from previous SDSS studies, while the choice of 
zo = 0.0725 minimizes the errors in k-corrections. 



3.2 Structural parameters 

The grizYJHK images were processed with 2DPHOT (La Bar- 
bera et al. 2008a) (hereafter LdC08), an automated software envi- 
ronment that allows several tasks, such as catalog extraction (us- 
ing S-Extractor), star/galaxy separation, and surface photometry to 
be performed. The images were processed using two Beowulf sys- 
tems. The optical images were processed at the INPE-LAC clus- 
ter facility, running 2DPHOT simultaneously on 40 CPUs. A num- 
ber of 31, 112 best-calibrated frames were processed in each of 
the griz wavebands, requiring ~ 2 days per band. The UKIDSS 
frames were processed at the Beowulf system available at INAF- 
OAC. A total of 12, 963 multiframes were processed by running 
2DPHOT on 32 CPUs, simultaneously. The processing took half a 
day for each band. 

A complete description of the 2DPHOT package can be found 
in LdC08; here we only outline the basic procedure followed to 
measure the relevant galaxy parameters. Both the optical and NIR 
images were processed with the same 2DPHOT setup to guar- 
antee a homogeneous derivation of structural parameters from g 
through K. For each frame, the so-called sure stars are identi- 
fied from the distribution of all the detected sources in the FWHM 
vs. signal-to-noise, S/N, diagram. This procedure allows an es- 
timate of the average seeing FWHM of the image to be obtained 
(see sec. 3 of LdC08). For each ETG, a local PSF model is con- 
structed by fitting the four closest stars to that galaxy with a sum 
of three two-dimensional Moffat functions. Deviations of the PSF 
from the circular shape are modeled by describing the isophotes 
of each Moffat function with Fourier-expanded ellipses. Galaxy 
images were fitted with PSF-convolved Sersic models having el- 
liptical isophotes plus a local background value. For each galaxy, 
the fit provides the following relevant parameters: the effective 
(half-light) radius, r c , the mean surface brightness within that ra- 
dius, < n >c,, the Sersic index (shape parameter) n, the axis ra- 
tio b/a, and the position angle of the major axis, PA. The total 
(apparent) magnitude, mr, of the model is given by the definition 
nix = — 2.51og(27r) — 51og(r e )+ < n> e - Mean surface bright- 
ness values are k-corrected and corrected for galactic extinction as 
described in Sec. 3.1. Moreover, cosmological dimming is removed 
by subtracting the term 10 log(l + z), where z is the SDSS spec- 
troscopic redshift. 

The characterization of the galaxy isophotal shape is done 
through the two-dimensional fitting of each ETG in the gri wave- 
bands, where Sersic models having isophotes described by Fourier- 
expanded ellipses are adopted. Only the fourth order cos term of 
the expansion, 04, is considered (boxiness - 04 < and disky- 
ness - 04 > 0, see e.g. Bender & Mollenhoff 1987, hereafter 
BM87). Only the gri band images are analyzed, since 04 is usu- 
ally measured in the optical wavebands. Since the models are PSF- 
convolved, the method above provides a global deconvolved esti- 
mate of (14. This estimate is somewhat different from the definition 
of BM87, where the peak value of 04 is derived in a given radial 
range, with the minimum radius being set to four times the see- 
ing FWHM and the maximum radius to twice the effective radius. 
Since many galaxies in the SPIDER sample have effective radii 
comparable to a few times the seeing FWHM value (see below), 
the BM87 procedure is not applicable. 

Fig. 4 shows the distribution of the average seeing FWHM 
value for all the retrieved images from g through K. The seeing 
FWHM was estimated from the sure star locus (see above). For 
each band, we estimate the median of the distribution of FWHM 
values, and the corresponding width values, using the bi-weight 
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Figure 4. Distribution of seeing FWHM values for all the grizYJHK 
frames (from left to right and top to bottom). The median value, fi, of each 
distribution is marked by the dashed line in each panel. Both the value of /i 
and the width, a, of the distributions are reported in the top-right corner of 
each plot. 



statistics (Beers, Flynn, & Gebhardt 1990). The median and width 
values are reported in Fig. 4 for each band. As expected, the me- 
dian FWHM value tends to smoothly decrease from the blue to 
NIR wavebands, varying from 1.24" in g-band to 0.82" in K-band. 
This variation corresponds to a relative change of ~ 34% (with re- 
spect to g band). Notice also that in YJHK (griz) bands almost 
all frames have seeing FWHM values better than 1.5" (1.8"), with 
90% of the values being smaller than 1.2" (1.5"). This decreasing 
of the seeing FWHM from g through K matches almost exactly the 
relative change of effective radii from optical to NIR wavebands 
(see below), making the ratio of FWHM to r e almost constant from 
g to K. 

Fig. 5 compares the distribution of x 2 values obtained from 
the two-dimensional fitting of galaxies in each band. The peak and 
width values of the distributions are computed by the bi-weight 
statistics and reported in the figure. The x 2 is computed as fol- 
lows. For each galaxy, we select only pixels la above the local sky 
background value. The intensity value of each pixel is computed 
from the two-dimensional seeing-convolved Sersic model. For the 
selected pixels, we compute the \ 2 as the rms of residuals between 
the galaxy image and the model. Residuals are normalized to the 
expected noise in each pixel, accounting for both background and 
photon noise. Notice that this \ 2 computation is somewhat differ- 
ent from that of the two-dimensional fitting procedure, where the 
sum of square residuals over all the galaxy stamp image is min- 
imized (see LdC08). This explains the fact that all the peak val- 
ues in Fig. 5 are slightly larger than one. An eye inspection of the 
residual maps, obtained by subtracting the models to the galaxy 
images, shows that the above \ 2 estimate is better correlated to 
the presence of faint morphological features (e.g. spiral arms, disk, 
etc.), that are not accounted for by the two-dimensional model. 
This is shown in Figs. 6 and 7. Both figures show residual maps in 
the r-band. Fig. 6 displays cases where the % 2 is c l° se to tne P ea k 
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Figure 5. Same as Fig. 4 for the \ 2 distributions. 



value (x 2 < 1.5), while Fig. 7 exhibits cases with higher \ 2 value 
(1.5 < \ 2 < 2.0). In most cases, as the \ 2 value increases, we can 
see some faint features to appear in the residual maps. We found 
that the percentage of galaxies with \ 2 > 1-5 is not negligible, 
amounting to ~ 17% in r-band. Most of the morphological fea- 
tures are expected to be caused by young stellar populations, hence 
disappearing when moving to NIR bands, where the galaxy light is 
dominated by the old, quiescent stars. From Fig. 5, one can actu- 
ally see that NIR bands exhibit a less pronounced tail of positive 
X 2 values with respect to the optical. This is also confirmed by a 
Kolmogorov-Smirnov (KS) test. For instance, in the case of r and 
K bands, the KS test gives a probability smaller than 1% for the 
corresponding x 2 distributions to be drawn from the same parent 
distribution. 



3.3 Uncertainty on structural parameters 

We estimate the uncertainties on structural parameters by compar- 
ing the differences in log r e , < /i > e , and log n between contiguous 
wavebands. To obtain independent estimates of the uncertainties 
on SDSS and UKIDSS parameters, the comparison is performed 
for the r and i bands, and the J and H bands, respectively. The 
variation of logr e , < /i > e , and logn with waveband depends 
on the measurement errors on structural parameters as well as on 
the intrinsic variation of stellar population properties (e.g. age and 
metallicity) across the galaxy, which implies a change of the light 
profile with waveband. As it is well known, this change is responsi- 
ble for the existence of radial color gradients inside ETGs (Peletier 
et al. 1990). The basic assumption here is that the r and i (J and 
H) bands are close enough that the variation of galaxy properties 
from one band to the other is dominated by the measurement er- 
rors. Tab. 1 shows that this is the case. In this table we report the 
sensitivity of color indices with respect to age and metallicity be- 
tween contiguous wavebands for a Simple Stellar Population (SSP) 
model. The sensitivities to age and metallicities, indicated as At 
and Az, are defined as the derivatives of the color indices with 
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Table 1. Sensitivity of color indices of an 
SSP model to age and metallicity. 



color index 
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0.313 
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Figure 6. Two-dimensional fit results for galaxies in r-band with typical 
X 2 value (x 2 < 1.5). Each plot shows the galaxy stamp (left) and the 
residual map (right) after model subtraction, using the same gray scale of 
intensity levels. The spatial scale is shown in the bottom-left corner of the 
left plots. For each galaxy image, the corresponding celestial coordinates 
(right ascension, RA, and declination, DEC, in degrees) and \ 2 value are 
reported in the upper-left corner. 




Figure 7. Same as Fig. 6 but for galaxies with high \ 2 value (1.5 < \ 2 < 
2.0 in r-band). 



respect to log t and log Z, where t and Z are the age and metal- 
licity, respectively. The derivatives are estimated as in La Barbera 
& de Carvalho (2009). We use an SSP from the Bruzual & Chariot 
(2003) synthesis code, with solar metallicity, Scalo IMF, and an age 
of t = 10.6 Gyr, corresponding to a formation redshift of Zf — 3 in 
the adopted cosmology. For the SDSS wavebands, the sensitivities 
reach a minimum in r — i, while they have a maximum in g — r, as 
expected by the fact that this color encompasses the 4000A break in 
the spectrum of ETGs at the median redshift of the SPIDER sample 
(z = 0.0725). The values of At and Az can be used to estimate 
the expected intrinsic waveband variation of structural parameters. 
As shown by Spolaor et al. (2009), bright ETGs have a large dis- 
persion in their radial metallicity gradients, Vz, with Vz varying 
in the range of about to —0.6 dex. On the contrary, age gradients 
play a minor role (La Barbera & de Carvalho 2009). Even consid- 
ering a dispersion of 0.3 dex in the metallicity gradients of ETGs, 
from Tab. 1, one can see that the corresponding scatter in the r — i 
internal color gradients would be only 0.079 ■ 0.3 ~ 0.024mag. 
Following Sparks & Jorgensen (1993), for a de Vaucouleurs profile, 
this implies an intrinsic scatter in the difference of r- and i-band ef- 
fective radii of only ~ 2.4%, hence much smaller than the typical 
measurement error on r e (see below). For the UKIDSS wavebands, 
the lowest sensitivities to age and metallicity are obtained in J—H, 
being even smaller than those of the optical colors. 

The errors on structural parameters are expected to be mainly 
driven by two parameters, the FWHM + pixel scale of the image 
(with respect to the galaxy size) and the signal-to-noise ratio, SN. 
Consequently, we bin the differences in log r e , < fi > e , and log n 
between r and i (J and H) bands with respect to the logarithm of 
the mean effective radius, logr e , and the SN per unit area of the 
galaxy image, SN/r 2 . The SN is defined as the mean value, be- 
tween the two bands, of the inverse of the uncertainties on Kron 
magnitudes. Each bin is chosen to have the same number of galax- 
ies. In a given bin, we estimate the measurement errors on logr e , 
< ii > e , and logn from the mean absolute deviation of the cor- 
responding differences in that bin. We refer to the measurement 
uncertainties as c logre , °"< M > e i an d f] 0g „, respectively. Since the 
errors on effective parameters are strongly correlated, we also de- 
rive the corresponding covariance terms, COVy,x, in each bin, 
where X and Y are two of the three quantities logr e , < pi > e , 
and log n. For each pair of X and Y, we perform a robust linear 
fit of the differences in Y as a function of the corresponding dif- 
ferences in X. The quantity COVy,x is obtained from the slope, 
s, of the fitted relation as COVy,x — s x a\, with the constraints 
\COV Y ,x K cry x a x , and COV Y ,x = COV x ,y- Fig. 8 plots the 



6 The robust regression is performed by minimizing the sum of absolute 
residuals of the Y vs. the X differences. 
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Figure 8. Errors on optical effective parameters as a function of the logarithm of the signal-to-noise (SAO per pixel. In the left panels, from top to bottom, we 
plot the uncertainty on the effective radius, the effective mean surface brightness, and the Sersic index. The left panels, from top to bottom, plot the covariance 
terms of the uncertainties on logr^, < fi > e , and logn. Different colors show different bins of the effective radius as shown in upper-right of the top-left 
panel. For a given color, the points shows the values of the uncertainties and the covariance terms obtained in different bins of the logarithm of the SN per 
pixel. The dashed curves are the best-fit functional forms used to model the dependence of the uncertainties on SN (see the text). 



quantities a. r , o" <)J>e , and <r log n as well as the relevant covari- 
ance terms as a function of log(SN/r^), for the r and i bands. The 
errors on effective parameters are strongly correlated to the SN per 
unit area. For log(S7V/r e )<2, as the signal-to-noise per unit area 
decreases, the errors tend to become larger. For \og(SN/r^)>2, 
all galaxies have effective radii comparable or even smaller than 
the pixel scale, with large values of log(SN/r'j) corresponding to 
smaller effective radii. As a result, the error on the effective param- 



eters tends to increase as well. The a, 



exhibits a similar behav- 



ior to that of a, and ov,,-, 

log r e <M> e 



though with a larger dispersion at 



given log(SN/ri). For the errors on effective parameters (<r k 



, , and COV. a 



) and Sersic index (a, ), the trends 



exhibited in Fig. 8 are well described by the following, empirical 
functional form: 



(2) 



,a^, a ndCOV iosrc ^ >c 



where y is one of the quantities cr lo r 
x — log(SA//re), P2 is a second order polynomial function, while 
ci and C2 are two parameters describing the increase of the error 
values at high SN per unit area. For the quantities COV <jl> x n 
and COV lo n lo , the trends in Fig. 8 can be modeled by a fourth 
order polynomial function. We have fitted the corresponding func- 
tional forms by minimizing the sum of absolute residuals in y. The 
best-fit curves are exhibited in Fig. 8. 

Fig. 9 plots the errors on structural parameters as derived with 
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Figure 9. Same as Fig. 8 but comparing J and H band structural parameters. 



the above procedure for the J and H bands. The trends are simi- 
lar to those obtained for the optical parameters, but with a larger 
dispersion, at given SN per unit area, which is likely explained 
by the fact that the number of galaxies in each bin with available 
photometry in J and H is smaller (by a factor of ~ 8) than that 
in r and i. The uncertainties on the NIR parameters are on average 
larger than those in the optical. For instance, at log(SN/r 2 ) — 2, 
the uncertainty on logr e is ~ 0.06 dex in the optical, increasing 
to ~ 0.1 dex in the NIR. The median errors on logr e , < fi > e , 
and logn, amount to ~ 0.1, ~ 0.5 mag/arcsec 2 , and ~ 0.1, re- 
spectively, in the optical, and to ~ 0.15, ~ 0.6 mag/arcsec 2 , and 
~ 0.12 in the NIR. This difference can be qualitatively explained 
by the fact that (i) the stars used for the PSF modeling have, on av- 
erage, a lower signal-to-noise ratio in the NIR than in the optical, 
and (ii) that the ratio of galaxy effective radii to the pixel scale of 
the images is smaller in the NIR than in the optical. The trends in 



Fig. 9 are modeled with the same functional forms as for the optical 
data. 



We use the above analysis to assign errors to the structural 
parameters to each galaxy in the SPIDER sample, for each wave- 
band. For a given galaxy, we first calculate its SN per unit area and 
then use the best-fitting functional forms to assign <r log rj , a , <ll>c , 
cr log n , and the corresponding covariance terms. For the gri bands, 
we adopt the functional forms obtained from the r — i comparison, 
while for the JHK bands we adopt the values obtained from the 
comparison of J- and ff-band parameters. In the z and Y bands, 
we apply both the optical and NIR functional forms, and then de- 
rive the errors by interpolating the two error estimates with respect 
to the effective wavelength of the passbands. 
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4 COLOR-MAGNITUDE RELATIONS 

As a first step in the comparison of optical to NIR properties of 
ETGs, we start to analyze the differences in their integrated proper- 
ties, i.e. the color indices. The goal is comparing total Sersic mag- 
nitudes in the different wavebands, in order to (i) characterize the 
completeness of the SPIDER sample in the space of the (Sersic) 
effective parameters (Sec. 8), and (ii) select suitable samples of 
ETGs for the analysis of the FP (see papers II and III). Thus, we 
estimate the color indices using the Sersic total magnitudes, rather 
than aperture magnitudes as in most of previous studies, and refer to 
them also as the total galaxy colors. This is by itself an important 
issue that will be addressed in a future contribution dealing with 
the different ways of measuring colors (de Carvalho et al. 2010, in 
preparation) which certainly goes beyond the scope of this paper. 

The comparison of color indices is performed by construct- 
ing the color-magnitude diagrams for different pairs of wavebands. 
Since the ETG sample has photometry available in eight wave- 
bands, we can derive seven different color-magnitude relations. We 
consider galaxy colors in the form of g — X with X — rizYJHK, 
and we write the color magnitude relations as: 

9-X = a gX +b gX X, (3) 

where a gX and b gX are the offsets and slopes of the relations. To 
simplify the notations, we also set a gg = and b gg = 0. Fig. 10 
plots the CM diagrams for all the 5, 080 ETGs with available data 
in all wavebands. In order to derive a v and b „ , we first bin each 
g — X vs. X diagram with respect to the magnitude X. We adopt 
N = 15 bins, with each bin including the same number of galax- 
ies. Varying the number of bins in the range of 15 to 25 changes the 
slope and offset values by less than 1%. For each bin, we derive the 
peak of the corresponding distribution of galaxy colors by applying 
the bi- weight location estimator (Beers, Flynn, & Gebhardt 1990). 
The uncertainty of the peak value is estimated as the standard de- 
viation of the peak values obtained in 1000 bootstrap iterations. 
This procedure has the advantage of being insensitive to outliers in 
the color distribution. The binning is performed up to a magnitude 
limit X, obtained by transforming the 2DPHOT r-band complete- 
ness limit of —20.55 (Sec. 8.1) through the median values of the 
g — X total colors. The values of a x and b x are then derived by 
fitting the binned values of g — X vs. X, with an ordinary least- 
squares fitting procedure with g — X as dependent variable. The 
uncertainties on a gX and b gX are obtained by randomly shifting 
(N — 1000 times) the binned values of g — X according to their 
uncertainties. The values of a gX and b gX are exhibited in Fig. 10. It 
is interesting to notice that, using total colors, the CM relations are 
essentially flat, with the slopes b gX being mostly consistent with 
zero within the corresponding uncertainties. In particular, the value 
of b x is consistent with zero at less than 2a in g — r, g — i, g — J, 
g — H, and g — K, while in g — z and g — Y the slopes differ from 
zero at ~ 3<j and ~ 5u levels, respectively. We do not find here 
any systematic steepening of the CM relation when enlarging the 
waveband baseline from g — r to g — K, as expected if the total CM 
relation would be purely driven by a mass-metallicity relation in 
ETGs. Following Scodeggio (2001), we can explain this surprising 
result by the fact that we use total color indices. ETGs have neg- 
ative color gradients, with color indices becoming bluer from the 
galaxy center to its periphery. As a result, when adopting colors in 
a fixed aperture, one is measuring the color inside a smaller region 
(with respect to r e ) for the brightest (hence larger) galaxies than for 
faintest galaxies in the sample. This leads to a misleading steepen- 
ing of the CM relation. We notice that using de Vaucouleurs model 



Table 2. Median peak values, a g -x, of the 
distributions of total colors, g — X (X = 
rizYJHK). The uncertainties are l<x stan- 
dard errors on the median values. 



color 


peak value 


9- 


r 


0.852 ± 0.007 


9 ~ 


- i 


1.274 ±0.004 


9 - 


- z 


1.589 ±0.008 


9 - 


Y 


2.282 ±0.012 


9 - 


J 


2.806 ±0.012 


9 ~ 


H 


3.472 ±0.012 


9 - 


K 


3.854 ±0.018 



magnitudes from SDSS would essentially lead to the same effect, 
as model magnitudes are estimated in a fixed aperture for all the 
SDSS wavebands (see Stoughton et al. 2002), and the SDSS effec- 
tive radii tend to be more underestimated (with respect to the Ser- 
sic r e ) for brighter than for fainter galaxies (see Fig. 16 of Sec. 6). 
Since the values of b gX — are mostly consistent with zero and our 
aim here is that of relating total magnitudes among different wave- 
bands (rather than performing a detailed study of CM relations), we 
have decided to set b gX — 0, and derive a g _ x as the median of the 
g — X peak values in the different magnitude bins. The values of 
a g X , together with the corresponding uncertainties, are listed in 
Tab. 2. The uncertainties are the errors on the median values. They 
are estimated from the width of the distribution of median values 
for 200 bootstrap iterations. 

Fixing the limiting magnitude of the ETG's sample in a given 
band, W (W = grizYJHK), one can use the CM relations 
to map that limit into equivalent magnitude limits, X Um (X = 
grizYJHK), in all wavebands. From Eq. 3, one obtains: 

a ... — a „ 1 + b ... 
X Um = -ttt-^ + ttt^ Wi im . (4) 
1 + b gX 1 + b gX 

In the particular case of b x — 0, one obtains the simplified ex- 
pressions, Xnm — a Y — a gX + Y, which is used in paper II to 
analyze the FP relation for color-selected samples of ETGs. 



5 DISTRIBUTION OF STRUCTURAL PARAMETERS 
FROM g THROUGH K 

Figures 1 1 to 13 exhibit the distributions of 2DPHOT Sersic param- 
eters from g through K. For each band, we select all the galaxies 
available in that band (see Sec. 2). Each distribution is characterized 
by its median value, //, and the width, a, estimated by the bi-weight 
statistics (Beers, Flynn, & Gebhardt 1990). Both values, fi and a, 
are reported in the plots. 

Fig. 1 1 compares the distributions of effective radii. The most 
noticeable feature is that the median value of logr e decreases 
smoothly from the optical to the NIR, varying from ~ 0.53 dex 
(r e ~ 3.4") in g to ~ 0.38 dex in K (r e ~ 2.4"). This change 
of ~ 0.15 dex corresponds to a relative variation of ~ 35% in 
r e , and is due to the fact that ETGs have negative internal color 
gradients, with the light profile becoming more concentrated in the 
center as one moves from shorter to longer wavelengths. Had we 
used the peak values of the distributions, estimated with the bi- 
weight statistics, rather than median values, the relative variation in 
r e would have been ~ 31% (~ 0.13 dex) instead of ~ 35%. The 
optical-NIR difference in log r e is in agreement with the value of 
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Figure 10. Color-magnitude relations for all the 5, 080 galaxies with photometry available in grizY J H K . From left to right and top to bottom, galaxy colors 
in the form of g — X, with Z = rizYJHK, are considered. For each panel: the slope and offset of the CM relation are reported in the upper-left corner, with 
the corresponding 1<t standard uncertainties; the colored circles mark the peak of the color distribution in magnitude bins, with the corresponding Icr error 
bars; the solid line shows the best-fitting CM relation; the dashed line mark the magnitude cut ion that band that corresponds to a magnitude limit of —20.55 
in the r-band (see text). 



- 0.125 dex (~ 29%) reported by Ko & Im (2005) for the sam- 
ple of 273 ETGs from Pahre (1999), between the V and K bands. 
Notice also that the relative change in r e fortuitously matches the 
improvement in average seeing FWHM between the g- and K-band 
images (Sec. 3.2), which makes the measurement of structural pa- 
rameters across the SPIDER wavelength baseline even more homo- 
geneous. 

Fig. 12 compares the distributions of values of the axis ra- 
tio, b/a, of the best-fitting Sersic models. The median as well as 
the width values of b/a turn out to be essentially constant from g 
through K, amounting to ~ 0.7 and ~ 0.2, respectively. The con- 
sistency of the b/a distributions with wavebands is in agreement 
with that found by Hyde & Bernardi (2009) when comparing g- 
and r- bands b/a values from SDSS. As expected, the fraction of 
ETGs decreases dramatically at low values of b/a with only a few 
percent of galaxies having axis ratios as low as 0.3. 

Fig. 13 compares the distributions of the Sersic index, n. Since 
we have selected bulge-dominated galaxies (fracDev r > 0.8, 
see Sec. 2), all the objects exhibit a Sersic index value larger than 
one, i.e. no galaxy has an exponential (disk-like) light profile. In 
particular, the fraction of ETGs becomes significantly larger than 
zero above the value of n ~ 2, which, according to Blanton et al. 
(2003b), roughly corresponds to the separation limit between blue 
and red galaxies in the SDSS. The distributions show a large scatter, 
with n ranging from ~ 2 to ~ 10. The median value of n is around 
6 for all wavebands, without any sharp wavelength dependence. 
On the other hand, we see some marginal change in the shape of the 
distribution with waveband. In the optical gri wavebands, a peak in 
the distribution is evident around n ~ 4. The distributions become 
essentially flat in the other (NIR) wavebands, with the exception 
of J where a peak is still present at n ~ 4. Some caution goes in 
interpreting these changes in shape. First, notice that observations 
and data reduction in the I-band of UKIDSS-LAS are carried out in 
a somewhat different manner with respect to the other wavebands 



-g ba 

U£ 


id 






I 

a 


i i i i | . 
= 0.53: 

=0.41 : 


-r bar 


'd 


_ 




i 

A 
c 


1 1 1 I - 

1=0.50: 

r=0.38 : 


-i ban 

L£ 


d 






1 

A 
c 


1 1 1 1 | . 

.=0.51 : 
=0.39 : 


-z bar 
.1 t-JT 


& 

d 






A" 
a 


=0.49: 
=0.43 : 

iItpi : 


-Y ba 
\l 




-t 
-+ 




^ 

C 


1=0.41 : 
r=o.36 : 

TV~P— t- 1 : 


-J bar 








A 
c 


.=0.41 : 
=0.34 : 


-H ba 

\I 




-t 

_L 


+- 


A 1 

G 
1 


=0.39: 

=0.37 : 


-K t 
", ["7 


a 






+- 


c 
I 


t=0.38: 
r=0.39 : 


1 2 

log r e (arcsec) 



12 12 

log r e (arcsec) log r e (arcsec) 



Figure 11. Distributions of logr e for all the grizYJHK wavebands 
(from left to right and top to bottom). The median value, /1, of each dis- 
tribution is marked by the vertical dashed line. The /j and width values 
(see the text) are reported in the upper-right corner of each panel. 

(see Warren et al.2007). A micro-stepping procedure, with integer 
pixel offsets between dithered exposures is performed. Images are 
interleaved to a subpixel grid and then stacked. This procedure re- 
sults in a better image resolution of 0.2" /pixel, with a better ac- 
curacy of the astrometric solution (useful for proper motion's mea- 
surements). We cannot exclude that this difference in data reduc- 
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Figure 12. Same as Fig. 1 1 for the distribution of axis ratios. 
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Figure 13. Same as Fig. 1 1 for the distribution of Sersic indices. 



tion affects the J-band distribution of n values. Moreover, as seen 
in Sec. 3.3, uncertainties on structural parameters change from g 
through K, hence preventing a straightforward comparison of the 
shape of the distributions among different wavebands. 

In order to analyze the waveband dependence of the 04 param- 
eter, we compare the 0,4 values among contiguous wavebands for 
the sample of 39,993 ETGs. The comparison is performed using 
only the gri wavebands, where the estimates are derived (see 
Sec. 3.2). Fig. 14 plots the differences in 054, 8(14, between g and r, 
and r and i, as a function of the median 04 value, < 04 >. The £04 
values are binned with respect to < 04 >, with each bin including 
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Figure 14. Differences of 04 estimates between r — i (upper panel) and 
g — r (lower panel) as a function of the median 04 value, < 04 >, in 
the gri wavebands. The blue dashed line marks the value of zero. The red 
curves are obtained by binning the data with respect to < 04 > . 



the same number (JV = 200) of galaxies. For a given bin, the me- 
dian difference of 8(14 values is computed. The median values are 
plotted as a continuous curve in Fig. 14, showing that there is no 
systematic trend of 0,4 from g through i. 



6 COMPARISON OF SDSS AND 2DPHOT STRUCTURAL 
PARAMETERS 

We compare the effective parameters measured with 2DPHOT with 
those derived from the SDSS photometric pipeline Photo, that fits 
galaxy images with two-dimensional seeing convolved de Vau- 
couleurs models (Stoughton et al. 2002). From now on, the dif- 
ferences are always in the sense of SDSS— 2DPHOT. The com- 
parison is done in r-band, by using the entire SPIDER sample 
of 39,993 ETGs (Sec. 2). Effective radii along the galaxy ma- 
jor axis are retrieved from the SDSS archive, and transformed 
to equivalent (circularized) effective radii, r e , with the axis ra- 
tio values listed in SDSS. The effective mean surface brightness, 

< /j, > e , is then computed from the circularized effective radii and 
the de Vaucouleurs model magnitude, m^ BV , r% using the definition 

< fi >e=mdev,r+5\og(r e ) + 2. 5 log(27r) . All magnitudes are 
dereddened for galactic extinction and k-corrected (see Sec. 3.1). 
For both 2DPHOT and Photo, we denote the effective parameters 
as r e and < fi > e , saying explicitly when we refer to either one or 
the other source. 

In order to compare the method itself to derive r e and < 
/i > e (rather than the kind of model, i.e. Sersic vs. de Vaucouleurs), 
we start by comparing the effective parameters of galaxies for 
which 2DPHOT gives a Sersic index of n ~ 4. To this effect, we se- 
lect all the ETGs with n in the range of 3.7 to 4.3, considering only 
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Figure IS. Differences in total magnitude, S mag (upper panel), effective 
radius, <5r e /r e (middle panel), and FP parameter, SFP (lower panel), be- 
tween SDSS and 2DPHOT, as a function of the SDSS r-band model mag- 
nitude. All values refer to the r-band. Red curves are obtained by bin- 
ning the differences with respect to m^ev r (see the text). The blue bars 
mark the range of expected differences due to the sky overestimation effect 
present in the SDSS parameters estimation. The size of the blue bars has 
been estimated from fig. 6 of Abazajian et al. (2009) for Smag and fig. 3 
of Adelman-McCarthy et al. (2008) for 8r e /r e . 



galaxies with better quality images (seeing FWHM < 1.5")- This 
selection results into a subsample of 4, 525 ETGs. Fig. 15 plots the 
differences in the total magnitude, effective radius, and FP param- 
eter as a function of mdev,r, where FP =logr- e — 0.3< n > e is 
the combination of logr e and < n > e entering the FP relation 
(Saglia et al. 2001). For the r e , we normalize the differences to the 
Photo r e values. For each quantity, the differences are binned with 
respect to m,iev,r, each bin including the same number (N = 200) 
of galaxies. In a given bin, the peak value of the distribution of 
differences (red curves in the figure) is computed by the bi-weight 
statistics (Beers, Flynn, & Gebhardt 1990). The SDSS total magni- 
tudes and effective radii differ systematically from those obtained 
with 2DPHOT, with total magnitudes (effective radii) being fainter 
(smaller) with respect to those of 2DPHOT. This effect tends to 
disappear for faint galaxies: the absolute differences in total mag- 
nitude decrease from ~ 0.15mag (~ 20%) at m<j CT ,r~ 15 to 
~ O.lmag (6%) at m,j e „ ir ~ 17.7. These differences are in the 
same sense as those reported from previous studies (e.g. Bernardi 
et al. 2007, Lauer et al. 2007), and can at least partly be accounted 
by the sky overestimation problem affecting SDSS model parame- 
ters (Adelman-McCarthy et al. 2008; Abazajian et al. 2009). Photo 
tends to overestimate the sky level near large bright galaxies, lead- 
ing to underestimate both total fluxes and effective radii. The effect 
has been quantified from the SDSS team by adding simulated r 1//4 
seeing-convolved models to SDSS images, and recovering their in- 
put parameters through Photo. Fig. 15 compares the range of val- 
ues for differences between input and output parameters from the 
SDSS simulations (blue bars), with what we find here. For r e , the 
average />/zoto-2DPHOT differences are only marginally consistent 
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Figure 16. Comparison of SDSS and 2DPHOT parameters as a function of 
the 2DPHOT Sersic index (left panels), and the SDSS total model magni- 
tude (right panels). From top to bottom, the same quantities as in Fig. 15 
are compared. The red curves are obtained by binning the data as in Fig. 15, 
but with each bin including 400 galaxies. 



with those expected from the simulations. For the total magnitudes, 
we find larger systematic differences, with 2DPHOT magnitudes 
being brighter, by a few tenths of mag, than what expected from 
the simulation's results. We should notice that larger differences in 
magnitude (and perhaps in r e ), in the same sense as we find here, 
have also been reported by D'Onofrio et al. (2008) when compar- 
ing their effective parameters with those from SDSS. Moreover, 
the SDSS simulations have been performed by assuming a given 
luminosity-size relation for ETGs, which might be slightly differ- 
ent for ETG's samples selected according to different criteria. Since 
one main goal of the SPIDER project is that of analyzing the FP re- 
lation, we have to point out that, although the differences in r e and 
rriT are significant, the FP quantity is in remarkable good agree- 
ment when comparing Photo to 2DPHOT. The average difference 
in the FP is less than a few percent and does not depend on the 
magnitude. As a consequence (see paper II), the FP coefficients 
change by only a few percent when using either 2DPHOT or Photo 
effective parameters. 

Fig. 16 compares differences between Photo and 2DPHOT 
parameters as a function of the Sersic index, n, as well as the 
(SDSS) absolute model magnitude in r band, °' 07 Mdev.r- Differ- 
ences are binned as in Fig. 15, considering only the TV = 39, 091 
galaxies with better image quality (see above). The comparison 
reveals large systematic differences, that strongly correlate with 
the Sersic index n. As n increases, Sersic total magnitudes be- 
come brighter - while Sersic effective radii become larger - than 
the SDSS values. The former trend is consistent with that reported 
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by Graham et al. (2005). Similar, but weaker, trends are also present 
as a function of the galaxy magnitude, when moving from fainter 
to brighter galaxies. This is somewhat expected, as ETGs exhibit a 
luminosity-Sersic index relation, with brighter galaxies having on 
average larger n (Caon, Capaccioli & D'Onofrio 1993). The trends 
of Fig. 16 are similar to those obtained by D'Onofrio et al. (2008) 
(see their fig. 1), when comparing Sersic effective parameters to the 
Photo quantities. As noticed above, the FP parameters are more 
stable with respect to the fitting procedure than the other quantities. 
In particular, the quantity FP shows only a weak dependence on 
galaxy magnitude (see lower-right panel of Fig. 16), with an end- 
to-end average variation of only 0.02 dex (~ 5%). Since the FP 
can be seen as a linear relation between the FP quantity and veloc- 
ity dispersion (or galaxy magnitude 7 , through the Faber- Jackson 
relation), the weak dependence of the FP quantity with magnitude 
implies that the coefficients of the FP are expected not to change 
significantly when using either 2DPHOT or SDSS parameters (see 
paper II). 



7 SPECTROSCOPY 

7.1 Velocity dispersions from SDSS and STARLIGHT 

We have re-computed central velocity dispersions for all the ETGs 
in the SPIDER sample. Velocity dispersions are usually measured 
by comparing the observed galaxy spectrum with single spectral 
templates, which are assumed to describe the dominant stellar pop- 
ulation of the galaxy. For ETGs, the spectra of red giant stars are 
usually adopted. On the other hand, ETGs frequently show mixed 
stellar populations, and this might significantly affect the cto es- 
timate for some fraction of the ETG's population. We derive the 
go's with the same procedure as in the SDSS-DR6 pipeline, i.e. 
the direct fitting of galaxy spectra (see Adelman-McCarthy et al. 
2008), but instead of using single spectral templates as in the SDSS 
pipeline, we construct a mixed-population spectral template for 
each galaxy. This is done as part of the automatic procedure de- 
scribed in STARLIGHT (Cid Fernandes et al. 2005) where velocity 
dispersion and stellar population parameters are determined simul- 
taneously. Hence, the new velocity dispersion values should be vir- 
tually unaffected by the different kinematics of the various stellar 
components. 

For each galaxy, we run the spectral fitting code STARLIGHT 
(Cid Fernandes et al. 2005) to find the combination of single stel- 
lar population (SSP) models that, normalized and broadened with 
a given sigma, best matches the observed spectrum (also normal- 
ized), which is first de-redshifted and corrected for extinction. We 
use SSP models from the MILES galaxy spectral library, with a 
Salpeter Initial Mass Function truncated at lower and upper cut- 
off mass values of 0.01 and 120 M , respectively (Vazdekis et 
al. 2010). These models are based on the MILES stellar library 
(Sanchez-Blazquez et al. 2006), which has an almost complete cov- 
erage of stellar atmospheric parameters at a relatively high and 
nearly constant spectral resolution of 2.3A (FWHM). This reso- 
lution is better than that of SDSS spectra, allowing us to suitably 
degrade the spectral models to match the resolution of the observed 
spectra (see below). 

One main issue for the estimate of velocity dispersions from 



7 Notice that the Sersic "n" is also correlated with velocity dispersion (Gra- 
ham 2002), but the correlation exhibits a large dispersion. 
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Figure 17. Wavelength dependent instrumental resolution for the SDSS 
spectra of five randomly selected ETGs. The resolution values are taken 
from the fits spectral files in the SDSS archive. The dashed line shows 
the resolution variation for the MILES SSP-models, which have a nearly 
constant value of 23 A (FWHM). For each galaxy, we degrade the MILES 
models (dashed curve) to match the corresponding wavelength dependent 
resolution (solid curves; see the text). 



SDSS is the wavelength variation of the SDSS spectral resolu- 
tion. For the ETG's spectra of the SPIDER sample, we found that 
the median value of the resolution varies from ~ 2.8A (FWHM) 
(o-inst ~ 90 km/s) in the blue (4000A) up to ~ 3.7A (FWHM) 
(o-inst ~ 55 km/s) in the red (8000A). Resolution also varies sig- 
nificantly among different spectra, as seen in Fig. 17, where we 
plot the a inat as a function of wavelength for the spectra of five 
randomly selected ETGs. Bernardi et al. (2003a) have accounted 
for the wavelength dependence of the SDSS spectral resolution by 
modeling it with a simple linear relation. On the contrary, in the 
present study, we do not perform any modeling of o~i ns t . For each 
galaxy in the SPIDER sample, we degrade the SSP models to match 
the wavelength-dependent resolution, ainst (\), of the correspond- 
ing spectrum. The a inst is measured from the SDSS pipeline by 
using a set of arc lamps, and provided in one fits extension of the 
spectrum fits file (see Stoughton et al. 2002). The MILES models 
are degraded by the transformation: 

M s (\) = J M{x)G{x-\)dc (5) 

where M(\) is a given SSP model, M s (\) is the smoothed model, 
and the function G(A) is a Gaussian kernel whose width is ob- 
tained by subtracting in quadrature the MILES resolution to the de- 
redshifted resolution Oi ns t (A/(l + 2)), where z is the galaxy spec- 
troscopic redshift. The integral is performed by discrete integra- 
tion. Notice that Eq. 5 reduces to a simple convolution in the case 
where a inst is a constant. For each galaxy, we run STARLIGHT 
using the corresponding smoothed MILES models. We use a set 
of SSP models covering a wide range of age and metallicity val- 
ues (with fixed solar [a/Fe]=0 abundance ratio). Age values range 
from 0.5 to ~ 18 Gyr, while metallicity values of log Z/Zq — 
-1.68, -1.28, -0.68, -0.38, 0., 0.2 are considered, resulting in a 
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Figure 18. Example of the advantage in using mixed-population templates 
to derive the galaxy central velocity dispersion. The upper panel shows an 
example ETG spectrum from SDSS (J154615.45-001025.4) and the best- 
matching synthetic models obtained by running STARLIGHT (i) with a set 
of 132 different SSP MILES models (red color), and (ii) with a single SSP 
model, having an age of 12.6 Gyr and solar metallicity (blue color). The 
(to values obtained in the two cases are reported in the upper panel. The 
lower panel shows the residuals obtained by subtracting the models to the 
observed spectrum. The gray bands in the upper panel show the masked 
windows we use to exclude from the fitting those regions possibly contam- 
inated by emission lines (see the arrows in the plot) as well as corrupted 
pixels (e.g. those corresponding to bad-columns). 

total of 132 SSP models. Fig. 18 illustrates the advantage of using 
a mixed-population rather than a single stellar population template. 
In order to exacerbates the difference between the two approaches, 
we selected the spectrum of one ETG for which STARLIGHT mea- 
sures a significant contribution from young (age< 2 Gyr) stellar 
populations. The Figure plots a portion of the spectrum, along with 
two best-fitting models obtained by either the mixed-population ap- 
proach (in red) or by running STARLIGHT with a single old stellar 
population template having an age of 12.6 Gyr and solar metallicity 
(in blue). Residuals are plotted for both cases. It is evident that the 
mixed-population model yields a better description of the contin- 
uum and the absorption features in the galaxy spectrum. In particu- 
lar, one may notice that the Mgb band (A ~ 5170A), which is one 
of the main spectral features to measure the <jq, shows significantly 
smaller residuals in the case of the mixed-population fit. In fact, the 
o"o value (reported in the upper panel of Fig. 18) changes dramat- 
ically from one case to the other. Fig. 18 also shows the masked 
regions used to avoid either corrupted spectral regions (e.g., bad 
columns) or regions possibly contaminated by nebular emission. In 
particular, we show three masked regions. The ones at A ~ 4850 
and 4980 A avoid the H/3 (4861 A) and [OIII] (4959 and 5007) 
emission lines, respectively, while the one at ~5220 A excludes 
pixels contaminated by a bad-column, as flagged in the SDSS spec- 
trum. 

Fig. 19 compares velocity dispersion values obtained from the 
SDSS spectroscopic pipeline with those measured in this work us- 
ing STARLIGHT. A good agreement is found, with only a small 
systematic trend at the low (< 90km/s) and high (> 280km/s) ends 
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Figure 19. Differences in logcro between STARLIGHT and SDSS as a 
function of the SDSS velocity dispersion value. Curves with different colors 
are obtained by binning the data with respect to cro(DR6), and taking for 
each bin the corresponding median value (red), the location value estimated 
through the the bi-weight statistics (magenta), and the 2.5cr clipped mean 
value (green). Notice that the trend with a is essentially the same for all 
the three estimators. The horizontal blue dashed line marks the value of 
zero, corresponding to a null difference between SDSS and STARLIGHT 
(jo values. 

of the (Jo range. In particular, for ao < 90km/s, the STARLIGHT 
velocity dispersions are slightly higher, by a few percent, with re- 
spect to those of the SDSS. It is important to emphasize that al- 
though STARLIGHT does not normalize the spectrum by the con- 
tinuum, which is done by the direct fitting method used in the SDSS 
pipeline, the good agreement found here is likely reflecting the ex- 
cellent quality of flux calibration obtained in DR6 (see Figure 7 
of Adelman-McCarthy et al. 2008). The impact of the above sys- 
tematic difference in ao on the scaling relations of ETGs is investi- 
gated in paper II. 



7.2 Uncertainties on the velocity dispersions 

To estimate the uncertainties on the STARLIGHT ao values, we 
looked for ETGs in the SPIDER sample having spectra with re- 
peated observations in SDSS. Out of all the 39,993 galaxies, we 
found 2,313 cases with a duplicate spectrum available. For all these 
duplicate spectra, we measure the ao with STARLIGHT, with the 
same setup and set of model templates as for the primary spec- 
tra. In each case, we compute the relative difference of ao as 
8a /a = (cto, 2 — o"o.i)/(To,i, where the indices 1 and 2 refer to the 
spectra with higher and lower S/N ratios. Fig. 20 plots the 8a /a as 
a function of the minimum S/N ratio, min(S/N), of each pair of 
duplicate spectra. The S/N is computed from the median S/N ratio 
in the spectral region of the Kg feature, within a window of 100A, 
centered at A = 4860 A. We bin the Sa/a values with respect to 
min(S/N), with each bin including 30 galaxies. For each bin, we 
compute the median and rms values of 5a /a. As the Fig. 20 shows, 
the median values (red color) are fully consistent with zero, imply- 
ing that, as expected (see Sec. 2), the S/N ratio of the spectra is large 
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Figure 20. Relative differences of ctq values, Sa/a, for ETG's spectra with 
repeated measurements available from SDSS. Each relative difference is 
computed with respect to the value obtained for the spectrum with highest 
S/N ratio. On the x-axis, the minimum S/N ratio, min(S/N), of each pair 
of spectra is plotted. The red solid curve shows the median difference values 
obtained by binning the Sa/cr's with respect to min(S/N), with each bin 
including the same number of 30 galaxies. The red dashed curve marks the 
value of zero in the plot. The blue solid curve shows rms of difference's 
values in each bin. The dashed blue line is obtained by modeling the rms 
values by Eq. 6. 



enough to obtain unbiased velocity dispersion estimates. The rms 
values, rmSa , provide an estimate of the error on STARLIGHT <to 
values. The rmscr increases at low S/N ratio and is well described, 
as shown in Fig. 20, by the following functional form: 



= 0.0383 + 3.2 x min(S/N)~ 



(6) 



For each galaxy in the SPIDER sample, we assign the error on 
the STARLIGHT ao value using the above equation, replacing 
min(S/N) with the median S/N ratio in the Hp region (see above) 
of the corresponding galaxy spectrum. Notice that this procedure 
assumes that the error on ao depends only on the S/N ratio, and 
is the same, for a given S/N, for both low and high ao galax- 
ies. In fact, we verified that considering only galaxies with ao < 
130 km s , the dashed blue curve in Fig. 20 does not change sig- 
nificantly, shifting upwards by less than 2%. Fig. 21 compares the 
distribution of uncertainties of the SDSS and STARLIGHT ao val- 
ues. The STARLIGHT distribution is significantly narrower than 
that of SDSS. This is likely due to the fact that we are using a 
one parameter (i.e. the S/N ratio) function to assign the errors (see 
Eq. 6). However, the peak position of both distributions is very sim- 
ilar. In fact, the median values of SDSS errors and STARLIGHT 
rmSa 's are fully consistent, amounting to ~ 0.07dex (~ 15%). 



8 COMPLETENESS 

8.1 SDSS vs. 2DPHOT completeness in r band 

We first analyze the completeness of the SPIDER sample in r band. 
The completeness of SDSS data is well characterized with respect 




Figure 21. Distribution of uncertainties on SDSS (dashed-black) and 
STARLIGHT (gray) ao values. A small offset of +0.004 has been applied 
to the gray histogram to make the plot more clear. 



to the Petrosian magnitude, m Ptr , which is one of the main selec- 
tion criteria used to target objects for SDSS spectroscopy. Here, 
we adopt m PiT = 17.77 as the reference value for the spectro- 
scopic completeness (Strauss et al. 2002). On the other hand, ef- 
fective parameters are defined in terms of either the model mag- 
nitudes (SDSS) or Sersic total magnitudes (2DPHOT). As shown 
in Blanton et al. (2001), the difference between Petrosian and 
model magnitudes depends on galaxy half-light radii. Moreover, 
as discussed in Sec. 6, there are large systematic differences be- 
tween Photo model magnitudes, rridev, and the 2DPHOT Sersic 
total magnitudes, tot. This implies that the m p limit maps into a 
different completeness magnitude value when using either SDSS 
or 2DPHOT parameters. In the following, we denote the absolute 

and total 2DPHOT magnitudes 



rridev magnitude as 07 M ( 



as M r , characterizing the completeness of the SPIDER sam- 
ple with respect to both om M d e V , r and 07 M r . We also use the 
symbol M r to refer indistinctly to either one or the other total mag- 
nitude. Fig. 22 plots SDSS and 2DPHOT absolute magnitudes as 
a function of m PtT . At a given m PlT , the scatter seen in the plot 
reflects all the different factors that enter the definition of total ab- 
solute magnitudes, i.e. the k- and galactic extinction corrections, 
the redshift range of the sample, as well as the intrinsic difference 
between M r and m p (see above). Assuming the SDSS sample to 
be complete down to m Pl r= 17.77, we define here 95% com- 
pleteness limits for Mdav,r and °' 07 M r , by adopting the ge- 
ometric approach illustrated in the figure. The method is similar 
to that described by Garilli, Maccagni, & Andreon (1999), where 
the completeness limit of a galaxy sample is defined as that magni- 
tude where galaxies begin to be missed in it because of the surface 
brightness detection limit. In practice, we select a range in absolute 
magnitude where all galaxies have m Pjr smaller than the SDSS 
spectroscopic completeness limit of 17.77. We consider the ranges 
of -20.7 to -21.7 and -21 to -22 for 07 M dOT , r and °' 07 M r , 
respectively. Then, in these ranges, we bin the distribution of m p , r 
values with respect to M r , and derive the 95% percentile of the 
m Pt r distribution in each bin. The binned points are fitted with 
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Table 3. Completeness limits in 
grizYJHK. Column 1: waveband. 
Column 2: 95% completeness magnitude. 
Column 3: number of galaxies brighter than 
the completeness limit. 
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consider here only 2DPHOT Sersic magnitudes, which are linearly 
related to logr e and < (x > e , and hence are the relevant quan- 
tities to characterize the completeness of the sample in the space 
of effective parameters. Fig. 23 plots the total Sersic magnitudes, 
°' 07 Mx, as a function of m p , r , where X is one of the available 
wavebands (i.e. X = grizYJHK). For a given band, we define 
an absolute magnitude range of °- 07 M Xl to °- 07 Mx 2 , where the 
sample is complete with respect to m v>r , and bin the corresponding 
distribution of m p>r values with respect to the °' 07 Mx ■ The 95% 
completeness magnitude in the band X is then defined as described 
in Sec. 8. 1 . For each band, the values of °- 07 M x x and °- 07 M x 2 are 
obtained by trasforming those in r-band through the median values 
of the ETG's color indices (see Tab. 2). For each band, the 95% 
completeness limit is reported in Tab. 3, together with the number 
of SPIDER ETGs brighter than that limit. 



m pr 

Figure 22. Completeness limits of the ETG sample for the SDSS absolute 
model magnitudes (upper-panel) and 2DPHOT absolute Sersic magnitudes 
(lower-panel) in r-band. The vertical dashed gray line in each panel marks 
the Petrosian magnitude limit, m p , r = 17.77, of the SDSS spectroscopy. 
The red curve is obtained by a second polynomial fit of the 95% percentile 
values of the m Pir distributions in different bins of total magnitude (see 
the text). The horizontal dashed lines mark the points where the red curves 
cross the vertical line, defining the 95% completeness limits of the sample. 



a second order polynomial, shown by the red curves in Fig. 22. 
The M r value where the polynomial intersects the vertical line 
of m p , r = 17.77 defines the point below which at least 95% of 
the galaxies, for whatever value of M r , are included in the sam- 
ple. We refer to the value of M r at the intersection point as the 
95% completeness limit of the sample. Notice that the 2DPHOT 
completeness limit is brighter than that of the model magnitudes. 
The difference, ~ 0.2mag, matches the average difference between 
' 0, ' Mdev,r and 0,07 M r measured for the faintest galaxies in the 
SPIDER sample (see Sec. 6). 

8.2 Completeness from g through K 

We apply the geometric approach of Sec. 8.1 to characterize the 
completeness limit of the SPIDER sample in all wavebands. Since 
our main goal is that of selecting volume limited samples for an- 
alyzing the FP relation in each waveband (papers II and III), we 



9 SUMMARY 

The data presented in this paper represents the most extensive ever 
obtained for early-type galaxies. Optical data from SDSS were en- 
tirely reprocessed; images with 2DPHOT (La Barbera et al. 2008) 
and spectra with STARLIGHT, as described above. This allows us 
to perform a proper comparison of different pipelines. We have cre- 
ated a database to store all of the photometric and spectroscopic 
parameters measured as part of our reprocessing. The SPIDER 
database (SdB) currently contains only the photometric (SDSS 
and UKIDSS) and spectroscopic (SDSS) data described here, but 
we expect to incorporate additional data from other wavelength 
regimes. The SdB architecture is extensible and designed to sup- 
port this growing process. We are currently developing an intuitive 
Graphic User Interface (GUI) which will allow the user to easily 
retrieve all information available in the SdB issuing SQL queries. 
Until the GUI is ready for general use we are making the data avail- 
able through the link 

http://www.lac.inpe.br/bravo/arquivos/SPIDER_data_paperI.ascii. 8 
The data presented here will be used in the forthcoming pa- 
pers analyzing the FP relations (papers II and III). Given the im- 
portance of obtaining meaningful measurements of the structural 
parameters entering the FP, this contribution examines how consis- 
tent these parameters are and their errors. Here we summarize the 
main properties of the sample defined in this paper and the charac- 
teristic parameters of each ETG: 

8 A mirror is also available at http://www.na.astro.it/ labarber/SPIDERA 
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Figure 23. Completeness limits in the grizY J H K bands. Each panel is the same as the lower panel in Fig. 22, but for a different waveband, X. The X runs 
from g to K going from left to right, and top to bottom in the figure. The upper-mid panel is the same as the lower panel of Fig. 22, and is repeated here for 
a better comparison among all the wavebands. The 95% completeness limits (see the text) are reported in the upper-left comers of the panels. The two solid 
gray horizontal lines indicate the absolute magnitude range of 07 Mxx to 0,07 Mx 2 , where the sample is complete with respect to m p>r , 



1 - When matching optical+NIR data we end-up with 5,080 
ETGs which can then be used to study the global properties of el- 
liptical galaxies in the nearby Universe (z<0. 1). 

2 - For each of these ETGs we have measured aperture and to- 
tal magnitudes, k-corrected and dereddened for galactic extinction. 
We also quantify for each galaxy the following essential parame- 
ters: effective (half-light) radius, r c , the mean surface brightness 
within this radius, < /i > Cj , the Sersic index (shape parameter) n, 
the axis ratio b/a, and the position angle of the major axis, PA. 
Besides, the galaxy isophotal shape is characterized by the a.4 pa- 
rameter. 

3 - Uncertainties in all of the parameters previously listed are 
estimated and presented as a function of the logarithm of the signal- 
to-noise per pixel. Median errors in r c , < fi > e ,, and n are ~ 0.1, 
~ 0.5 mag/arcsec 2 , and ~ 0.1, respectively, in the optical, and 
~ 0.15, ~ 0.6 mag/arcsec 2 , and ~ 0.12 in the NIR. 

4 - We do not find any systematic steepening of the CM rela- 
tion when considering a waveband baseline from g — r to g — K, 



as expected if the total CM relation is solely driven by a mass- 
metallicity relation in ETGs. According to Scodeggio (2001), this 
result can be interpreted as follows: ETGs have negative color gra- 
dients, with color indices becoming bluer from the galaxy center 
to its outskirts. Thus, when we adopt colors in a fixed aperture, we 
measure a color inside a smaller region (with respect to r e ) for the 
brightest (hence larger) galaxies than for faintest galaxies in the 
sample, leading to a misleading steepening of the CM relation 

5 - The study presented here of the structural parameters of 
ETGs reveals the following properties: log r c decreases by 35% 
from the optical to the NIR, reflecting the internal color gradients 
in these systems; the axis ratios, b/a, have their distribution essen- 
tially constant from g through K, with a media value of ~0.7 and 
width of ~0.2; The Sersic index is always larger than 1, since we 
selected only bulge dominated galaxies and spans a domain from 
~2 to ~10, with a median of 6 for all wavebands; no systematic 
trend of a4 was found from g to i. 

6 - We present measurements of central velocity dispersion 
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obtained using STARLIGHT. Extensive comparison with the esti- 
mates provided by the SDSS pipeline, shows that our estimates are 
unbiased over the whole a range, with a median error of ~15%. 

7 - Comparison of two independent pipelines was done (SDSS 
and 2DPHOT) and reveals significant differences in magnitude and 
effective radius with respect to the Sersic index and absolute model 
magnitude in r-band. The impact of such differences on the scaling 
relations of ETGs will be addressed in paper II. 
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